WChPT analysis of twisted mass lattice data 
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We perform a Wilson Chiral Perturbation Theory (WChPT) analysis of quenched twisted mass 
lattice data. The data were generated by two independent groups with three different choices for 
the critical mass. For one choice, the so-called pion mass definition, one observes a strong curvature 
for small quark masses in various mesonic observables ("bending phenomenon"). Performing a 
combined fit to the next-to-leading (NLO) expressions, we find that WChPT describes the data 
very well and the fits provide very reasonable values for the low-energy parameters. 
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I. INTRODUCTION 



Twisted mass lattice QCD |l|,|2j has many advantages 
for numerical lattice simulations, with automatic 0(a) 
improvement at maximal twist H 0, being proba- 
bly the most striking one (for a recent review see Ref. 
[6j). Maximal twist is achieved by tuning the bare un- 
twisted mass mo to a particular (critical) value such that 
some matrix element vanishes. The condition for max- 
imal twist is not unique and various choice have been 
employed in quenched simulations 0, H, U El El • 

A puzzle observed in early quenched simulations is the 
so-called "bending phenomenon" 0] ■ Employing the pion 
mass definition for maximal twist (i.e. tuning mo to the 
value where the pion mass would vanish without a twisted 
mass term) , one observed a strong curvature in the quark 
mass dependence in many observables (m^, f^,m p ) for 
small twisted quark masses /i. This unexpected obser- 
vation spurred further numerical simulations with other 
definitions for maximal twist 0, H E3 as well as theoret- 
ical studies. Nowadays the bending phenomenon seems 
well understood both in terms of the Symanzik effective 
theory E2 as well as m Wilson Chiral Perturbation the- 
ory (WChPT) [H [lljll (for introductions to lattice 
ChPT see also Refs. ^MMfl). 

It is a pleasant side effect of this effort to understand 
the bending phenomenon that there are lots of data avail- 
able for five different lattice spacings and three definitions 
of maximal twist. Moreover, light quark masses could be 
reached with values for /m p down to ~ 0.3, where one 
might expect WChPT to provide an effective description 
of the data. In particular, the characteristic curvature of 
the bending phenomenon provides a distinctive test for 
WChPT to pass if it is the correct low energy effective 
theory for twisted mass lattice QCD. Provided WChPT 
describes the data very well we may also obtain estimates 
for some combinations of low energy constants of the ef- 
fective theory. This was sufficient motivation for us to 



carry out a WChPT analysis of the existing lattice data. 
Preliminary results involving data for two definitions of 
maximal twist only can be found in Ref. jl8| . 



II. FITTING THE DATA 

We analyzed quenched lattice data generated by two 
different groups 

II H IS 113 with the Wilson plaque- 
tte action at B = 5.85 (a » O.I23fm) and (3 = 6.0 
(a « 0.093fm). The standard Wilson fermion action with 
a twisted mass term was employed to calculate a variety 
of mesonic observables. The twisted quark mass cov- 
ers the range m w « 270 ... 1200 MeV, or, alternatively, 
m^/mp ss 0.3 — 0.8. Besides the pion mass there is data 
available for the pseudoscalar decay constant, the angle 
cott^wT, as well as some more observables which we did 
not analyze. The untwisted bare quark mass was tuned 
according to three different definitions of maximal twist: 
the pion mass definition, the PCAC mass definition and 
the parity conservation definition. In total there are at 
most 52 data points available for each lattice spacing. 

There exists a vast literature on WChPT for twisted 
mass lattice QCD El El El HI El, which contains ex- 
pressions for various mesonic observables up to next-to- 
lcading order (NLO). The lattice artifacts are included 
through order 0(a 2 ) for different power countings and 
definitions of maximal twist. Here we use the NLO for- 
mulae given in Ref. EHj which include all NLO terms 
consistently for the two regimes where either it ~ a or 
fi ~ a 2 . The quenched chiral logarithm [22, |23| is also 
included in the formulae. 

We performed combined fits of the WChPT formulae 
at NLO to the data for three observables: 
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where V® and A a ^ denote the (nonsinglet) vector and ax- 
ial vector current, respectively. At NLO we have in total 
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FIG. 1: Results of a combined fit for f v and R at /3 = 6.0 (xVd.o.f. = 0.23 with d.o.f. = 26). The upper plots show the results 
where all data points are included in the fit, while data points with afxo < 0.0302 only are included in the fits shown in the 
lower plots. 



thirteen free fit parameters. Even though this is a fairly 
large number it is still small compared to the number of 
data points. 

We performed various fits, starting with all data points 
included and then successively remove the data points at 
high quark masses. In all cases we obtain good fit results 
with x 2 /d.o.f. w 0.2 — 0.5, even if all data points up to 
m^/trip w 0.8 are included. 1 Since we do not trust ChPT 
to work at such high masses we prefer to drop the data 
for the highest three masses. The fit results for and 
R at (3 = 6.0 are shown in figure ^ Even in this fit the 
heaviest point corresponds to m^/mp w 0.63, which is 
still heavy. Dropping more data points, however, makes 
the fit more and more unstable, so we cannot reduce the 
number of data points much further. 

Apparently WChPT describes the data very well. In 
particular the bending for small masses in case of the pion 
mass definition is very well reproduced. This feature is 
independent of the number of data points included in 
the fit, even though the values for the fit parameters are 



different (see below). Note that the curvature in the data 
for R with a/io > 0.3 is also well described even though 
the heavier data points are excluded from the fit. 

The fits give reasonable values for the fit parameters. 
For the quenched chiral log parameter Sq, for example, 
we find 
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0.10±0.03 
0.054 ±0.011 



(3 = 6.0 
/3 = 5.85 



(3) 



which is in very good agreement with the results obtained 
by other groups (for a summary see Ref. |24|). 

We also obtain an estimate for the low-energy constant 
c-2- This parameter was first introduced in Ref. [25| and 
enters the chiral Lagrangian according to 2 



L x = 



(4) 



The sign of C2 determines the phase diagram of the lattice 
theory and the pion mass splitting Am 2 = rn^ a — m^ ± 
in the chiral limit is given by C2// 2 |2*5j . 



1 Note that the \ 2 value is underestimated since the data is highly 
correlated. 



2 The definition of C2 is not uni que and it is sometimes defined 
differently, for example in Ref. iRfl . 
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From our fits we obtain for C2 the value 

c 2 = [291MeV+f^] 4 8 = 5.85 (5) 

This is the first determination of this low-energy con- 
stant, so we cannot compare with other results. However, 
the value seems reasonable based on dimensional analysis 
arguments. The fit for the smaller lattice spacing with 
(3 = 6.0 does not determine c 2 very well; for the mean 
value we obtain approximately [170 MeV] 4 , but the error 
is of the same size. 

The physical parameters, on the other hand, are very 
well determined by the fit. For the pseudoscalar decay 
constant in the chiral limit we find 

f 141.2McV± 1% 8 = 6.0 

fo = < (6) 
( 141.4MeV± 1% 8 = 5.85 

which is in very good agreement with earlier determina- 
tions. For the low-energy constant |26| . entering the 
NLO expression for the decay constant, we obtain 

f 1.03(5) 3 = 6.0 

«b = \ (7) 
[ 0.97(4) 8 = 5.85 

Also these values agree very well with previous results in 
Ref. j2|| . Note that the results for f and a§ do not show 
any significant dependence on the lattice spacing. This 
is expected if WChPT works, since the main dependence 
on a is captured by other terms in the chiral expansion, 
being directly proportional to powers in a. 

We emphasize that the errors we quoted so far are only 
the statistical errors given by MINUIT which we used to 
perform the fits. These errors are underestimated due to 
the highly correlated data and the true statistical error 



can be substantially larger. A second error source are 
systematic uncertainties, induced, for example, by the 
number of data points included in the fit. It is not simple 
to give a precise estimate for this error but we observed 
that the central value for /o changes by roughly 3 percent 
for the different fits we performed, while varies by 
about 7 percent. 

III. CONCLUSIONS 

We performed fits of WChPT to quenched twisted 
mass data for , /„. and the Ward-Takahashi angle 
cotwwT- We find that the NLO expressions describe 
the data very well with small \ 2 values and reasonable 
values for the low-energy parameters. In particular, the 
bending phenomenon in case of the pion mass definition 
is very well reproduced. 

The bending phenomenon is a very characteristic fea- 
ture of twisted mass lattice QCD. It is encouraging that 
WChPT describes this distinct curvature very well. This 
indicates that WChPT, i.e. ChPT for lattice QCD, seems 
to work. Previous studies [23, IH, |2(| , using untwisted 
Wilson fermions, came to contradicting results and were 
not conclusive at all. 

So far we performed separate fits for each lattice spac- 
ing. In a next step it would be very interesting to per- 
form a combined fit to the entire data set and take the 
continuum limit. These results should be compared to 
the results one obtains after a standard continuum ex- 
trapolation where one assumes a polynomial lattice spac- 
ing dependence. This would partly answer the question 
whether WChPT is not only able to describe the lattice 
data but also necessary to extract the correct continuum 
physics from the data. 
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